Spin Correlation and Discrete Symmetry in Spinor Bose-Einstein Condensates 
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We study spin correlations in Bose-Einstein condensates of spin 1 bosons with scatterings dom- 
inated by a total spin equal 2 channel. We show that the low energy spin dynamics in the system 
can be mapped into an o(n) nonlinear sigma model(NL(TM). n = 3 at the zero magnetic field limit 
and n = 2 in the presence of weak magnetic fields. In an ordered phase, the ground state has a 
discrete Z2 symmetry and the degeneracy space is [(7(1) x S"~^]/Z2. We explore consequences of 
the discrete symmetry and propose some measurements to probe it. 



Recently, there has been a burst of theoretical and 
experimental activities on spinor Bose-Einstein conden- 
sates An optical trap confines alkali atoms inde- 
pendent of spins and liberates spin degree of freedoms 
For sodium(^^A^a) or rubidium(*^_R6) atoms with 
nuclear spin 7 = 3/2 and electrons at s orbits, the energy 
splitting between hyperfine multiplets is of order lOOrnfc. 
At temperatures as low as lOOnk, '^^Na and ^'^ Rh atoms 
can be considered as simple bosons with a hyperfine spin 
F — 1. The ground state of N spin 1 noninteracting 
bosons has [N +\)N /2 folds spin degeneracy, by contrast 
to its magnetically trapped cousins. Presumably, hyper- 
fine spin-dependent scattering lifts the degeneracy and 
leads to a spin correlated state. Optically trapped BEG 
therefore sets up a platform for studying quantum mag- 
netism in many boson systems and adds a new dimension 
to already extremely rich physics in these systems. 

The spin-dependent two -body interaction in BEG is 
characterized by t/2(i"i,r2) = 5{vi — r2)[co + C2F1F2], as 
suggested in an early paper |Q. Here cq = (50 + "^92) 1^, 
C2 — (<?2 — 5o)/3; gp — ^TrTi^ap /M , M is the mass of the 
atom and of is the s-wave scattering length in the total 
spin F channel. Thus, the spin correlation in a BEG is 
determined by C2. For ^'^ Kb, 32 < 5o or C2 < and the 
scattering is dominated by the total spin F = channel. 
In the ground state, all spins of atoms prefer to align in a 
certain direction and have a maxima magnetization Q. 

For ^^JVa studied experimentally the scattering 
between '^^Na atoms is dominated by the total spin 
F ~ 2 channel, i.e. 32 > go- The scattering between 
'^^Na atoms thus leads to an " antiferromagnetic" spin 
correlation. Efforts have been made to understand the 
ground state properties, exact excitation spectra, collec- 
tive modes Many interesting predictions, such as 
spin waves, spin mixing dynamics were made for ^^Na 
BEG where C2 > 0. 

In this paper, we show the spin dynamics in BEG with 



C2 > is characterized by an 0(71) nonlinear sigma model 
(NLaM) of n components, n = 3 at zero magnetic field 
limit and n = 2 in the presence of a weak magnetic field. 
Spin correlations in spinor BEG can be studied in the 
context of the NLcrM. We identify that the internal or- 
der parameter space for BEG as [S^ x S'^]/2'2(zero field 
limit) and the ground state is nematically ordered. We 
explore consequences of the discrete Z2 symmetry. 

To describe the spin correlated BEG, it is most conve- 
nient to introduce Weyl representation of SU{2) involv- 
ing polynomials of a unit vector {u,v) Each unit 
vector is represented by a point SI on a sphere with po- 
lar coordinates [0,4)); namely u — exp{i(f>/2) cos{0/2), 
V — exp{~i(f)/2) sm6/2. The corresponding hyper- 
fine spin operators are F^ — ud/dv, F~ = vd/du, 
and Fz — {ud/du — vd/dv)/2. The scalar prod- 
uct between two wavefunctions g and / is defined as 
/ g*{u,v)f{u,v)(Kl/A-K. (We reserve Vl for the spin ro- 
tations discussed below.) 

Under spin rotations TZ ~ exp(z_FjXi/2) exp(zFy0i/2) 
exp(z_F'2(/)i/2), u and v transform into 

M(f^i,Xi) = exp(zxi/2) 

(cos 6*1/2 exp(—j(/)i/2)w + sin 0i/2 exp(i(/)i/2)w), 
v{^2,X2) = cxp(-iX2/2) 

(- sin 02/2exp(-i02/2)u + cos 612/2 exp(i02/2)w), (1) 

where fJi 2 — (^1,2, 0i,2)- Spin-1 wavefunctions are poly- 
nomials of degree 2 in u and v. %/3u^, \/&uv, \/?>v^ cor- 
respond to m = 1,0,-1 states. All F — \ states can 
also be expressed in term of \/Qu{Q,x)v{Q,-2) with r2i,2 
properly chosen. 

The Hamiltonian for spin-1 bosons can be written as 
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The second term is hyperfine spin-dependent interaction 
with C2 > and the last term is the couphng with an 
external magnetic field H = TJe^; 5 is a Lande factor of 
an atom and /xs is the Bohr magneton. 

The wavefunction of N spin-1 Bosons generally can be 
written as |^ 

^'({r„}) = 7'n„=i...jv$JV„(ra)V6Wa(r2la(ra))Wa(f22a(ra)). 

(3) 

V is the permutation of {N^}, {Oiq, r22a}- labels 
an one-particle orbital state. Phases xi,2 in Eq.l are 
absorbed by a gauge transformation of the wavefunction 
introduced above. For BEC under consideration, we 
take r2iQ,,2Q — ^^i,2(r) and $Ar^ (r) = $(r) ($(r) is a com- 
plex scalar field). By introducing n(r) — (fJi + r22)/2, 
L(r) = (r2i-r22)/2 and <I>(r) = yj p{y) exp(«x(r)), we de- 
rive from Eqs.2,3 an effective Hamiltonian as a function 
of two pairs of variables : {n(r),L(r)} and {p(r),x(r)}- 
These are collective variables of the N interacting spin- 
1 Bosons, which describe the spin dynamics and phase 
dynamics respectively. In this representation, the hyper- 
fine spin dependent interaction in Eq.2 is mapped into 
Hs = J (irc2L^(r)yO^(r), which only depends on collective 
variable L when p(r) is taken as a constant. This indi- 
cates that each atom acquires an inertial Iq = l/2c2/3 in 
the presence of hyperfine spin dependent scatterings with 
C2 > 0; the rotation energy in the presence of a finite spin 
moment L is thus h^/2Io. 

In the most interesting limit, we can introduce the lo- 
cal spin density as l(r) = L(r)p(r). n and 1 satisfy the 
constrain n(r) • l(r) = 0. Commutation relations between 
p and X: n(r) and l(r) are given as [yo(r), x(r')] = ihS{r — 
r'); [n„(r),n0(r')] = 0, [l„(r), n^(r')] = ine"l'''n^S{r ~ 
r'), [la(r),l,3(r')] = ihe°'''^n-y6{r - r'). e"^'^ is an anti- 
symmetric tensor. These identities are valid when L per 
atom is much less than unity and n(r) can be considered 
as a classical " vector" , components of which commutate 
with each other. The corresponding Lagrangian density 
can be derived a,s C = Cs + Cc + Csc with 
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/:, = -^[(Vn(r)f + l(9.nn, 
2m vf 

^cs^Y7^[nxdrn-nx{Wx-^)n]. (4) 

Here p ^ p(0), = ^Ipc^jM, Vs = y^2pc2/M. We 
introduce t ^ it as the imaginary time. Nonlinearality 
is imposed via a constraint |n^| = 1 at a low frequency 
limit. Eq.4 is the main result of the mapping and we have 
kept contributions which are lowest order in terms of dr 
and V. Cc is taken in a Gaussian approximation and 
should be replaced by a full Gross-Pitaveskii Lagrangian 
in general. We will be mostly interested in the spin sector 



and simplification in Cc doesn't affect conclusions here. 
Cs in Eq.4 represents an o(3) NLcrM. 

The last term Csc characterizes a coupling between a 
spin rotation and the superflow in BEC due to Berry's 
phase. This term, however, is linear in dr and quadratic 
in spatial gradient and is negligible compared with Cs , Cc 
at the long wave length limit. Particularly, such a cou- 
pling vanishes in a configuration where L is zero. 

At the zero field limit, there exists a saddle point solu- 
tion for the spin sector n(r) — no, l(r) — 0. ng lives on a 
unit sphere. By expanding Eq.4 around the saddle point 
solution, we obtain spin waves with sound like spectrum 
uj = y^4c2p/Mk, which can also be obtained in the Gross- 
Pitaveskii approach Q|. However, Eq.4 here is valid for 
any point fii on the unit sphere. The effective 

NLctM derived here allows us to describe spin correlated 
states well beyond the Gross-Pitaveskii approach. 

For a symmetry broken state, the ground state wave- 
function in Eq.3 with fJi = = n is invariant under a 
global transformation n, % — )■ — n, x -I- tt; 

^{n, x) = *(-n, X + tt), *(n) = (-l)^*(-n) (5) 

where x is the phase of the scalar field ^(x) intro- 
duced before. In obtaining this symmetry, we notice 
u{n) = exp(i7r/2)z;(— n), with 7r/2 from a phase of a spin- 
1/2 particle under a 180" rotation. Eq.5 shows that the 
many-body wavefunctions characterized by (n, x) and 

n,x + 7r) are indistinguishable. Thus, the spinor BEC 
under considerations has a discrete Z2 symmetry and the 
order parameter space is TZ = [S^ x S'^]/Z2, with Z2 as a 
two-element group of integers modulo 2. 

As in a classical uniaxial nematic liquid crystal where 
diatomic molecules are indistinguishable upon an inver- 
sion of their directors n and the internal order param- 
eter space is S'^/Z2, the Z2 symmetry identified here 
also indicates that there exists a tensor order parameter 
Q . For the purpose of demonstration, let us introduce 
a "director" d^; = [uv* + vu*]/2, dy — [uv* — vu*]/2i, 
dz = [uu* — vv*]/2. One then can show for the ground 
state wavefunction in Eq. 3, < d >= but nematic order 
parameter Qap = l/p[< d^d^ > — l/3Tr < d^d^ >] is 
nonzero. Here <> stands for an average taken over the 
ground state wave function In fact, 
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Qaf3 = ^[-3noQno/3 + Sap]. (6) 

According to Eq.6, the director d aligns in a plane per- 
pendicular to no. 

Topological defects in this case are of particular inter- 
est because of the discrete symmetry in Eq.5. Following 
the general theory for the classification of defects in a 
symmetry broken state, the Z2 symmetry leads to tt spin 
disclinations superimposed with half vortices, which oth- 
erwise don't exist. The corresponding wavefunction of 
composite linear singularities {Z2 strings) is 
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Here ^ = x+iy and lines are located at = xo+iyo', n, m 
are integers and is superfluid velocity [0 . Each string 
is characterized by {m,n). However, (— l,n), (±3,n), 
(±5, n) strings can be obtained by deforming string (1, n) 
and are homotopical identical to (1, n). 

In a composite string given in Eq.7, n evolves into — n 
along a loop enclosing ^o- The corresponding spin wave- 
function changes its sign under an inversion n ^ — n, 
following the identity u{ja)v{ja) = — m(— n)w(— n). A su- 
perflow Vs of a half vortex is present to compensate the tt 
phase under n ^ — n rotation and ensure the single val- 
uedness of the wavefunction. This composite structure is 
different from a linear defect in a classical nematic liquid 
where tt disclinations are free topological excitations. In 
fact, in a coherent spinor BEC, a bare 7r-spin disclination 
carries a cut along which phase changes abruptly from tt 
to 2ti or 0. This cut starting at the disclination ends 
only at the boundary of the system and costs an energy 
linear in term of the size of BEC. For a similar reason, 
the energy cost to have a tt disclination and a half-vortex 
separated at a distance L is linearly proportional to L. 
Composite strings in Eq.7 should be considered as results 
of confinement of tt spin disclinations and half vortices in 
spinor BEC |ll|. 

In the presence of an external magnetic field along z 
direction, Jli 2 = n ± e^g^^BH / Ac2P and n satisfies con- 
strain 
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Obviously, an external magnetic field breaks the S*^ sym- 
metry and confines low frequency sector of n in a plane 
perpendicular to H itself. The Lagrangian in the pres- 
ence of a magnetic field is that of an 0(2) NLctM; it has 
S^-symmetry at the frequency lo <^ psH. At the high 
frequency limit, n precesses in a field 4c2l(a;), much larger 
than the external field and S'^-symmetry is restored. 

As a consequence, the order parameter space for the 
quantum spin nematic state in an external magnetic field 
is 7?. = [S^ X S^]/Z2- The nematic order parameter Qap 
is still given by Eq.7, with the easy plane of d parallel to 
the external field. Wavefunctions for linear defects are of 
same forms as those in Eq.7, but with all strings (m,7i) 
homotopically distinguishable. 

We have restricted ourselves to the weak magnetic field 
limit and neglected the possible quadratic Zeemann shift 
Ti-QZ — X^a Q^^^za (the external field H is along di- 
rection) . Inclusion of the quadratic Zeemann shift yields 
an additional term Cqz = J pQH^{n'^ + n^)dx to the 



NLcrM derived. The main effect of this contribution is to 
align n along the external field. When this shift domi- 
nates, the ground state is left with a double degeneracy: 
n = and n = — e^. The spin wave develops an energy 
gap of order QH^. We will focus on the zero magnetic 
field case in the rest of discussions. 

In general, following Eq.4, spin correlated BEC can be 
studied by considering a NLctM, 
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We introduce dimensionless length and time: f = 
rp^/"^, T = TVsP^/^. Derivatives are defined as 
{df,dx,dy, di). is a square root of the ratio between 
potential energy at an interatomic scale p^^^ /2m and 
zero point kinetic (rotation) energy C2p/2 of an individual 
atom. The zero point motion is absent for noninteracting 
bosons but gets stronger when C2 increases, or the effec- 
tive inertial Iq gets smaller. The o(3) NLcrM has order 
and disordered phases at d > 1, depending on the param- 
eter /. Most of qualitative results about spin correlations 
in BEC can be obtained in a renormalization group (RG) 
approach ||]. For pAa^ = 10~^ as is in experiments, a 
long range nematic order should be observed. 

In ID case, at zero temperature and zero field, one 
should expect there will be no long range order and the 
state is nematically disordered, following the RG results 
of NLcrM jl^. These nematic disordered states mimic 
the quantum spin liquid states proposed in the literature 
of Heisenberg antiferromagnetic systems(HAFS). How- 
ever, we notice that the NLcrM derived from the micro- 
scopic Hamiltonian in this paper doesn't have a O-term 
Co = 0/Att J drdxii ■ dn/dr x dn/dx, which is gener- 
ally present in HAFS studied before jl^. Absence of 
a 0-term, which ensures an energy gap in the excitation 
spectrum of nematic state, follows a fact that the Berry's 
phase under rotations of n vanishes identically ^ . 

Two remarks are in order. 1) At a high density limit, 
one should also take into three-body, four-body elastic 
scatterings. This can further modify the short distance 
dynamics but will not affect conclusions arrived above in 
a qualitative way. 2)0ne should be cautious about the 
definition of " phase" since the alkali atoms under investi- 
gation are in a long lived metal stable gaseous state. The 
life time of alkali atomic gas is limited by three-body in- 
elastic collisions . The collision rate is proportional to 
the square of the number density of atoms and increases 
dramatically as the density increases. This sets a prac- 
tical limit in order for the quantum disordered nematic 
liquid to be probed in BEC. 

Since the experiment |Q was done in a BEC cloud with 
a few millions sodium atoms, it is also particularly inter- 
esting to consider the symmetry restoring due to a fi- 
nite size effect. We take a weakly interacting limit where 
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quantum fluctuations of finite wave length are negligible. 
In a zero mode approximation, Hamiltonian becomes, 



"Hz.m — pC2 



2N' 



(10) 



L is the total spin of the N-Bosons system and can be con- 
sidered as an angular momentum operator defined on a 
unit sphere of n p|. And [n, n] = 0, [L^, n^] = ihe°'^''n-y, 
[La,L^] — ihe'^^'L.y. The energy spectrum is given by 
L/^ = l{l + 1), which is identical to that obtained in 
I = 0,2, 4, ....N, if N is an even number; and Z = 1, 3, 5, ... 
otherwise. The energy gap between lowest lying excita- 
tions is inversely proportional of N and vanishes at large 
N hmit. 

Now consider a wave packet with n confined within 
a region centered at no — Bz of size V < S^n >o ^ 1 
on the unit sphere at t — 0. Spread < S^n > is de- 
fined as the expectation value of + ^ly- I^i spherical 
polar coordinates {9, (j)) where n is a vector represented 
by (sin 9 cos (j), sin 9 sin (j), cos 9) , the corresponding wave 
packet (for an even N) can be constructed as 
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(11) 



Here B ~ exp(— /^/2ct); < 5^n >o= Ao/a with con- 
stant Aq estimated to be a constant of unity. Yifl(9,(l)) 
are spheric harmonics, a ^ 1 but a/N is vanish- 
ingly small. The energy of this wave packet is AE — 
2AQac2p/N . In the limit N — s- oo, an infinitesimal ex- 
ternal field will stablize this wave packet with respect to 
the rotation invariant state. 
The spread of n at a time t is 



< {Snf >t=< (Snf >o +4:Ao(7{t 
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(12) 



which is valid at i < Tc = (t/AqAE. At t > Tc, the 
spread is of order unity. Therefore, < S^n >t exceeds 
the initial spread V < S'^n >o at a time of order 1/AE. 
At a longer time Tc, n starts to explore the whole unit 
sphere S*^ and the rotation symmetry is restored due to 
spin-dependent scatterings in BEC. 

Eq.l2 also imposes restrictions on measurements. A 
measurement of n in ^^iVa BEC discussed here excites 
the ground state to an excited state of energy AE, where 
n has a finite spread on and AE is infinitesimal small 
in the thermal dynamical limit. The time scale to re- 
store the broken symmetry is determined by the two- 
body spin dependent scattering lengths and the number 
of atoms in BEC. A measurement taken at a time scale 
longer than a/AoAE should reveal the symmetry restor- 
ing because of zero point motions of n. With 10^ 
and C2P ^ 100Hz{500nK), the symmetry restoring time 



is about one day which is much longer than the life time 
of the trap. However, for a smaller cloud with 10^"'* 
atoms, the symmetry restoring can take place within 
10 — lOOsecs, before recombination processes take place. 

Nematically ordered BEC has very fascinating optical 
properties |^. In the presence of spin-orbital couphngs, 
dielectric constant is a tensor expressed in terms of Qa[3- 
This suggests that birefringence takes place in the system 
as a direct evidence of the hidden Z2 symmetry. Another 
experiment consequence associated with the broken sym- 
metry is the enhanced small angle light scattering due to 
thermal fiuctuations of n. For '^^Na, the light scattering 
amplitude can be four order of magnitude higher than 
that in an isotropic BEC and the nematic BEC is opti- 
cally turbid. Finally, the Z2 symmetry also implies that 
there exists a local Z2 gauge field in BEC. This was re- 
cently considered and will be published elsewhere |]l5| . 
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